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presently available to accurately quantify the homogeneity of 
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computer and analyzed for various area fractions. 
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have been analyzed as well. A statistically sufficient number 
of particles for use in the model was determined and used for 
all simulations. Average first nearest neighbor spacing 
values for dilute arrays of particles approach those of 
Poisson distributions of infinitesimal points. As the 
particle density increases, the average spacing values 
approach those of hexagonal arrays. For low area fractions 
there is little distinction between random and non-random 
distributions, both from statistical and visual perspectives. 
For higher area fractions there is a discernable difference 
between the statistical data for random and non-random 
distributions, but the visual differences are more obvious. 
These observations hold for both constant size particles and 
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I . INTRODUCTION 

Most engineering structural applications involve the use 
of materials containing second phases in the form of discrete 
precipitates or distributions of particulates. This latter 
category includes discontinuously reinforced metal-matrix 
composites. The behavior of these materials during processing 
and the final mechanical properties depend upon the nature, 
volume fraction, size and also the dispersion of the 
particles . 

In both conventional alloys and in metal-matrix composites 
the second phase particles may be obstacles to the movement of 
dislocations. If a dislocation is prevented from completely 
traversing the crystal under an applied stress, then the 
crystal will become harder to deform. A higher stress is then 
required to move the dislocation past the obstacle and 
therefore the material must be stronger. Orowan [Ref. 1] 
found that the shear stress required for a dislocation to 
bypass a pair of small particles may be expressed as: 




where is the shear modulus of the matrix, b is the burgers 
vector and D is the space between particles in the slip plane. 
From this relationship it can be seen that a smaller inter- 
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particle spacing will result in a larger stress to move the 
dislocation past the particles. 

Larger second phase particles will also affect nucleation 
and grain growth during recrystallization of conventional 
alloys. In general it has been found that greater inter- 
particle spacing and larger particles will facilitate the 
process of recrystallization [Ref. 2]. Close particle spacing 
and small particles have the opposite effect. The second- 
phase particles increase the dislocation density when the 
material is plastically deformed, increasing the driving force 
for recrystallization. However, they also impede the 
rearrangement of dislocations to form mobil high angle grain 
boundaries and thereby hinder the migration of these grain 
boundaries, thus retarding recrystallization. 

In research on A1 alloys with an AlCUj particle size 
greater than 0.5 /xm it was found that inter-particle spacing 
had the greatest effect on recrystallization [Ref. 2]. For 
wide particle spacings, subgrains formed by the rearrangement 
of dislocations and grew until they reached a critical size, 
at which point they would begin to migrate. The movement of 
the subgrains was impeded by the particles. This allowed for 
further nucleation of additional subgrains and resulted in 
accelerated recrystallization. 

The deformation structure of a deformed alloy is also 
dependent on the particle size. Large particles cause local 
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distortion of the lattice and increase the rate of nucleation. 
Small particles increase the homogeneity of the dislocation 
distribution which slows down the nucleation rate. Therefore, 
for an arrangement of coarse, widely dispersed particles the 
recrystallization is significantly increased. On the other 
hand, for a small particle diameter and small inter-particle 
spacing, the dislocation distribution is much more dense and 
homogeneous which reduces the rate of recrystallization. 

Second phase particle distributions may be either random 
or non-random. Random or uniform distribution infers that 
local particle densities are the same throughout the material. 
Non-random distributions are just the opposite and may be the 
result of insufficient processing of the material. High and 
low density bands of particles may be the result of processing 
an MMC where the original material consisted of clustered 
groups of particles [Ref. 3]. 

Kocks [Ref. 4] showed that the average area a swept out by 
a single dislocation segment is dependent on the applied 
stress a. This relationship is shown, in Figure 1. It can be 
seen that there exist a stress in which a dislocation can 
sweep out an area and keep moving indefinitely. This 
asymptotic stress is the macroscopic flow stress a. It is 
also apparent that the greatest flow stress for a given 
particle density is for a regular array of particles, is least 
for a clustered array, and the flow stress for a random array 
lies between the two. Thus the distribution of second phase 
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particles has a significant role in determining the mechanical 
properties of two phase materials. 



o 




Figure 1. Applied stress a vs mean free slip area a of 
one dislocation loop, for strong obstacles of area density 
l/£^, in square, random and clustered arrays [Ref. 3]. 

Clearly, it would be desirable to have a means of 
accurately assessing the different types of microstructures 
resulting from various processing techniques. With this 
knowledge these techniques can be optimized to efficiently 
achieve the desired properties of the material. Current 
methods used to determine these values assume that the 
distribution is random and that the particles are 
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infinitesimal points, which is certainly not the case for real 
materials. Methods do exist for finite sized particles but 
here the assumption is that they are distributed in a regular 
geometric array. 

This study, which is essentially follows on from that of 
M. Pas [Ref. 5], was undertaken to quantify such values as 
average first nearest neighbor particle spacings (5) for 
various area fractions, particle size distributions and for 
non-random cases as well. Any trends or patterns in the 
associated data was investigated as well. In order to 
facilitate mathematical and computer analysis computations 
were kept on a two dimensional basis. 
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II . BACKGROUND 



The nearest neighbor separation distance between second- 
phase particles governs many properties of two phase 
materials. Properties such as strength, ductility, and 
toughness are sensitive to the mean separation distance and 
may also be affected by the range and uniformity of this 
distance throughout the material. Current methods for 
determining this spacing for particles of finite size is 
actually based on the assumption that the particles are 
infinitesimal points which do not occupy space. Throughout 
this thesis, unless otherwise specified, the word "point” 
shall be regarded as an infinitesimally small entity which 
does not occupy space. The word "particle" shall be defined 
as a small two or three dimensional object. 

Underwood [Ref. 6] has described a method for determining 
the average first nearest neighbor distance, S, for randomly 
distributed points. The method is based on the probability of 
finding a nearest neighbor to a given point within an annulus, 
of radius r and width dr, centered on the point (Figure. 2) . 
This probability can be expressed as: 



P(r) dr 



1 - rP{r) dr 

J 0 



2nNj^dr 



( 2 . 1 ) 



where N,i^ is the number of particles per unit area and r is the 
distance between a point and its first nearest neighbor. The 
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Figure 2. The probability of locating a nearest neighbor 
to a point may be found by integrating over an annulus 
centered on the point. 



function P(r) which satisfies equation (2.1) was shown to be: 

P(r) = (2.2) 

Substituting equation (2.2) into the expression that defines 
the average first nearest neighbor distance results in the 
following expression for S for a random distribution of 
infinitesimal points: 

6 = f“rP(r)dr = 0.500^’;^^^ (2,3) 

J 0 
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On the other hand, if particles or points are arranged in 
a regular hexagonal array, Figure 3, the nearest neighbor 
distance for each particle will be the same everywhere and 
will be equal to <Sh„ for the array. A geometrical analysis 
results in the following expression for 

6^^^. = 1.075^;"/" (2.4) 




Figure 3. Particles arranged in a hexagonal array. 

Based on equations (2.3) and (2.4) it can be seen that, 
for a given value of r “^Hex. is greater than <S for a random 
distribution. 

If particles (or points) are arranged into any regular 
array the first nearest neighbor distance (NND) , r, would be 
the same for each particle and therefore would equal the 
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resulting value of S for the array. If any one of these 
particles (or points) were shifted out of its position in the 
array it would come closer to at least one other particle. 
Now, there would be at least two particles with a lower first 
NND, but none with a larger first NND since the remaining 
particles would still have other neighbors at a distance of r. 
Therefore, as the particles are shifted out of a perfect 
geometric array and become more random the value of S will 
decrease towards a lower limit defined by the Poisson 
distribution. 

Distributions of finite size particles can not truly be 
random. The area (or volume) occupied by a particle removes 
the possibility that another particle can be positioned in a 
portion of that area. This is the effect of non-overlapping 
particles. However, for dilute arrays (low area fractions) 
the possibility of two particles overlapping is small. 
Therefore, it is reasonable to expect that the value of S for 
dilute arrays of finite size particles would in fact approach 
that of the Poisson distribution. 

In previous work, performed by Mike Pas [Ref. 5], attempts 
to demonstrate that dilute arrays approached the Poisson 
distribution were not completely successful. Pas plotted 
vs S for various area fractions of constant size particles. 
The plots asymptotically approached values of S greater than 
those of the Poisson distribution with decreasing area 
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fraction. However, Pas did show that as the area fraction of 
particles increases, values for 6 shift towards values 
approaching that of a hexagonal array, equation (2.4). Pas 
also found that in non-random (high and low density bands of 
particles) distributions, as the percentage of particles in 
the high density bands increases, the overall value for S 
decreases . 

The ultimate goal of this work would be an analytical 
solution to quantify the degree of randomness of the particle 
distributions in microstructures. Incorporating this solution 
with a computer based image analyzing system would make it 
possible to readily assess the mechanical properties dependent 
on the particle distributions. 
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III. EXPERIMENTAL PROCEDURE 

A. OVERVIEW 

The simulations of particle distributions were generated 
using a computer. Four independent programs were written, one 
for each of four cases analyzed: random distribution of 

constant-size particles (RANPARTl) ; non-random distribution of 
constant size particles (RANPART2) ; random distribution of 
particles with a log-normal size distribution (RANPART3) ; and 
non-random distribution of particles with a log normal size 
distribution (RANPART4) . RANPARTl, 2, 3 and 4 are included in 
appendices A, B, C and D respectively. 

For the first two cases the algorithms were based on 
I programs originally developed by M. Pas [Ref. 5] and written 
in TURBO PASCAL. Here, Pas's programs were rewritten in 
FORTRAN 77 and incorporated several other modifications as 
deemed necessary. The programs were run either on an IBM- 
compatible PC or a Digital Equipment Corp. VAX 3100 work 
station. 

B. RANDOM DISTRIBUTION OF CONSTANT SIZE PARTICLES, RANPARTl 
The flow chart for RANPARTl is shown in Figure 4. In the 

input portion of the program the user enters the following 
parameters: 
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SEGREGATE PARTICLES 
ALONG THE BORDER 
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DETERMINE NEAREST 
NEIGHBOR SPACING 
VALUES 




> 








OUTPUT DATA 





Figure 4. Flow chart for the program RANPARTlJOR 
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• Area percentage of particles for the defined region; 

• X dimension of the defined region; 

• y dimension of the defined region; 

• Minimum spacing between particles (if particles are 
allowed to come in contact, enter 0) ; 

• Enter "n" if a specified number of particles is required 
I or enter "r" if a specified particles radius is required; 

I • If entered "n", enter the total number of particles, if 
I entered "r", enter the particles radius; 

• Enter "A" if particles are of finite size or enter "p" if 
the particles are infinitesimal points. 

From these parameters the total area of the defined region, 

! total area occupied by the particles, the area per particle, 

and the total number of particles or particle radius 

(whichever value was not input) are calculated. 

The particles are positioned in the next portion of the 

: program. A pseudorandom number generator is called twice to 

provide x and y coordinates for a given particle. This 

particle is then checked to ensure that it does not overlap 

any other previously sited particle. This is accomplished by 

calculating the distance to every other particle. If this 

distance is less than twice the sum of the particle radius and 

the minimum spacing, the position is rejected and the program 

loops back to generate a new set of coordinates. This process 

is repeated until the total number of particles required to 

establish the desired area fraction has been met. 

Once all the particle locations have been established, 

particles located along the border are segregated from those 
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in the interior of the defined region. This is done so that 
errors in first nearest neighbor distance (NND) calculations 
associated with the border will be minimized. The underlying 
rationale for this step in the program will be discussed in 
the next chapter. All particle positions are checked and 
those having coordinates placing them within one particle 
radius of the border are moved to a separate array in the 
program memory. 

First NND calculations are made by calculating for each 
interior particle the distance to every other particle, 
including those located along the border, and then finding and 
saving the minimum first NND value. This process is repeated 
for every other interior particle. The average first NND, 5, 
is calculated by summing the minimum first NND values and 
dividing by total number of interior particles . 

In order to plot the distributions of first NND values an 
increment had to be determined. To do this, the Sturgis rule 
was applied as it is an established method for determining the 
number of class intervals from a set of values. The Sturgis 
rule is defined by the following expression: 

N = 1 + 3.31og(n) (3.1) 

Where N is the number of class intervals and n is the 
number of values in the set. The value of n for this work is 
equal to the number of interior particles. The number of 
class intervals was divided into the difference between the 
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maximum and minimum first NND values. This result (defined as 
"dpinc" in the programs) was used as the incremental range for 
sorting first NND values into equal sized intervals. The 
number of NND values in each interval was then divided by the 
total number of interior particles to get a relative quantity. 
The output of the program consists of all input parameters, 
the calculated values described above and a calculated value 
for Na, the number of particles per unit area. 

RANPARTl was run for area fractions ranging from 1% to 50% 
utilizing a particle radius of 0.5 units. The defined area 
was adjusted in each run to maintain a nominal value Np«5000 
particles. The minimum spacing was kept at 0 for all 
simulations . 

Each run was repeated using the same parameters, but now 
declaring the particles to be points. The results consisted 
of point distribution and a distribution of finite sized 
particles with the same value of N^. 

For each simulation equation (2.1) was used to plot a 
Poisson distribution corresponding the simulated distribution 
of particles and points. The value of N,^ came from the 
simulated point distribution. The probability density 
function was plotted as a function of r in the range from the 
minimum and maximum values of first NND derived from the point 
distribution. The calculated increment Ar, based on the 
Sturgis rule, was substituted into the equation for dr. 
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C. NON-RANDOM DISTRIBUTION OF CONSTANT SIZE PARTICLES, 

RANPART2 

The following modifications were made to RANPARTl to 
produce RANPART2 . The defined area in which the particles are 
positioned was divided into five sub-areas of equal size. In 
the input portion of the program the user declares what 
proportion of the total area occupied by the particles will 
be assigned to each subsection. The sum of the percentages 
must equal 100%. 

After each particle (for finite sized particles) is 
checked for overlap, the program determines which sub-area the 
particle (or point) belongs. By maintaining a running sum on 
the area occupied by the particles in each sub-area, the 
computer determines whether the addition of another particle 
in a given sub-area would cause the area occupied by the 
particles to exceed the limit. The particle position is 
rejected if the limit will be exceeded. If this happens the 
program loops back to generate another pair of particle 
coordinates. When the process is complete the desired degree 
of banding has been achieved. Figure 5 is the flow chart for 
RANPART2 . 

All runs made for RANPART2 used 10%, 35%, 10%, 35%, and 
10% sub-area particle densities. This was considered to 
provide a relatively large deviation between the high and low 
density bands. The simulations were made for area fractions 
ranging from 0.01 to 0.30. 
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Figure 6. Flow chart for the program RANPART2.FOR 
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D. RANDOM DISTRIBUTION OF PARTICLES WITH A LOGNORMAL SIZE 

DISTRIBUTION, RANPART3 

In RANPART3 a second random number generator is used to 
determine a particle radius from a lognormal distribution. 
The call for the particle radius is made before the call for 
the particle coordinates so that, if particle overlap occurs, 
only a new position will be generated but the existing radius 
will still be used. Only when a particle is finally 
positioned or, if particle overlap occurs 600 times 
consecutively for the same particle radius, is a new particle 
radius generated. This process is required in order to 
maintain a lognormal size distribution and to prevent biasing 
towards smaller sized particles as the particle density 
increases. The program must now use individual particle radii 
in subsequent calculations. 

The input portion of the program was modified as follows. 
The user must enter a value approximately equal to the natural 
logarithm of the mean particle radius desired, which is the 
mean of the underlying normal distribution, and a standard 
deviation. Several iterations of these two variables may be 
required to get the desired mean radius, especially for larger 
standard deviations. Since a point distribution derived from 
the same input parameters of constant size particles would be 
no different than that for particles of varying size, this 
option is not available for the programs with a lognormal 
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particle size distribution. The total number of particles 
entered in the input portion of the program must be of 
sufficient quantity to ensure that the number of interior 
particles generated is enough to meet the desired area 
fraction.) Since now the particles are no longer of constant 
size there is no direct relationship between the number of 
particles and the total area occupied by the particles. A 
running sum of calculated particle areas is maintained. 
Particle generation continues until this sum is equal to the 
value corresponding to the desired area fraction. 

Segregation of particles along the border is performed 
much the same; but here the largest particle radius in the 
distribution is used to establish the extent of the border. 
Any particle whose center lies outside this border is 
segregated . 

The maximum, minimum and mean particle radius are 
determined and provided at the output. The flow chart for 
RANPART3 is shown in Figure 6 . 

RANPART3 was run for the same area fractions as RANPARTl. 
Values of -3.07 and 0.4 were used as input for the mean of the 
underlying normal distribution and for the standard deviation, 
respectively. 
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Figure 6. Flow chart for the program RANPART3.FOR 
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E. NON-RANDOM DISTRIBUTION OF PARTICLES WITH A LOGNORMAL SIZE 

DISTRIBUTION, RANPART4 

RANPART4 was derived from RANPART3 in the same way 
RANPART2 was derived from RANPARTl. Five sub-areas were 
incorporated into RANPART4 just as they were incorporated into 
RANPART2. To maintain a record of the area occupied by the 
particles in each sub-area, the area of each new particle 
needs to be calculated and added to the existing area in which 
it lies. 

The flow chart for RANPART4 is shown in Figure 7. 
RANPART4 was run for the same area fractions as RANPART2 and 
for the same values of mean and standard deviation as 
RANPART3 . 




Figure 7. Flow chart for the program RANPART4 J*XDR 
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IV. 



RESULTS AND DISCUSSION 



I A. OVERVIEW 

The first phase of this v/ork consisted of a study to 
' determine a statistically sufficient number of particles to be 
used in the simulations. The problem of particles positioned 
I along the border of the array was also addressed in this 
phase. In the second phase, four general cases were studied: 

I random distribution of points and finite sized particles of 
constant size; non-random distribution of points and 
particles of constant size; random distribution of particles 
with a log normal size distribution; and finally, non-random 
I distribution of particles with a log normal size distribution. 

B. DETERMINATION OF A SUFFICIENT NUMBER OF PARTICLES AND THE 
PROBLEM OF PARTICLES POSITIONED ALONG THE BORDER 
In order to determine a statistically sufficient number of 
> particles to be used in the simulations the program RANPARTl 
was modified. The modifications essentially consisted of 
enclosing RANPARTl inside two loops. This allowed repeated 
i program running for the same parameters. Then, the average of 
I the results of repeated runs, was calculated. This ejqseriment 
was used to calculate an average value of 6 , over a given 
number of runs, for various numbers of particles (Np) , 

■ holding constant. The test was conducted for two area 
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fractions, 0.10 and 0.40 and for point distributions as well. 
Na was arbitrarily set equal to 100. The results of the tests 
are shown in Table 1 and in Figure 8. 



Np 


• 

O 

11 


af = 0.4 


Points 


No. 1 

Runs 1 


10 


0.06801 


0.08902 


0.05781 


1000 


20 


0.06562 


0.08633 


0.05556 


1000 


30 


0.06442 


0.08524 


0.05425 


1000 


40 


0.06376 


0.08455 


0.05362 


1000 


50 


0.06349 


0.08434 


0.05326 


1000 


60 


0.06349 


0.08400 


0.05282 


500 


70 


0.06282 


0.08370 


0.05271 


500 


60 


0.06281 


0.08352 


0.05257 


500 


60 


0.06281 


0.06349 


0.05240 


500 


<00 


0.06233 


0.08326 


0.05240 


100 


200 


0.06173 


0.08255 


0.05146 


100 


400 


0.06154 


0.08223 


0.05131 


100 


400 


0.0612T 


0.08211 


0.05106 


100 


500 


0.06127 


0.08198 


0.05107 


50 


600 


0.06107 


0.08195 


0.05083 


50 


700 


0.06099 


0.08191 


0.05098 


25 


800 


0.06091 


0.08172 


0.05046 


25 


900 


0.06087 


0.08172 


0.05087 


25 


1000 


0.06076 


0. 0617S 


0.05075 


5 


2000 


0.06099 


0.0819C 


0.05025 


5 


3000 


0.06035 


0.08434 


0.05059 


5 


4000 


0.06048 


0.08434 


0.05035 


5 


5000 


0.06041 


0.08127 


0.05030 


5 



Table 1. Results from the modified RANPARTl for area 
fractions of O.l and 0.4 and for points, = 100. 
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AVERAGE FIRST NEAREST NEIGHBOR DISTANCE 



0.10 



0.09 



0.08 



"“Theoretical 

•-Particles, a^ = 0.40, rp = 0.05 
4-Particles, a^ = 0.10, rp = 0.05 
^—Calculated Points 




0.07 — 



0.06 — 




♦ » ♦ 



0.05 




10 100 1000 10000 
TOTAL NUMBER OF PARTICLES 

Figure 8. Total number of particles vs average first nearest 
neighbor distance for “ 100 
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For the random distribution of points, S approaches the 
theoretical value of 0.05 for Np > lOOO. For both area 
fractions, 0.40 and 0.10, S approaches and appears to remain 
constant at values of about 0.081 and 0.060 respectively, 
again at Np > 1000 in both cases. 

In developing the method for positioning circles or points 
on a plane the dimensions of the planar surface had to be 
established. Within this finite sized area all the particles 
were positioned. From this situation it became evident that 
when calculating first nearest neighbor distances, a particle 
whose center lies on or near the border of the defined region 
could conceivably have a closer neighbor than that calculated. 
Figure 9 shows a defined area and a particle, a, located near 
the border. Particle a' s nearest neighbor within the area is 
particle b. However, it is possible that if the defined area 
had been larger, particle a may have a nearest neighbor even 
closer than b such as particle c. Clearly this would be the 
situation in any real material (imagine the field of view of 
a micrograph being expanded) . Since the existence of these 
additional particles could only serve to lower the first NND, 
any calculations made for 6 without taking this into 
considerations would be overestimated. 

Calculations for S were repeated, holding constant and 
for the same area fractions, as described above. The model 
was modified such that first NND calculations were not made 
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Figure 9. The border problem. Particlee positioned along the border of 
a defined area (a) may have a potential nearest neighbor (c), if the 
defined region was expanded, which is closer than the nearest neighbor 
positioned inside the area (b). 
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for particles whose centers fell within one particle radius 
from the border. These border particles were used in 
determining first NND for particles located inside this 
exclusion range. The results are shown in Figure 10. 

This attempt to eliminate the border problem was only 
partially successful as shown in Figure 10 and has the 
greatest effect for low numbers of particles. The values for 
S were consistently lower for calculations involving both 
particles and points when the numbers of points or particles 
utilized was below 1000. For Np greater than 1000 the effects 
of using a larger number of particles or points compensates 
for the border problem by simply increasing the relative 
number of particles located inside the defined area. 

The data compiled in Table 1 was used to determine a more 
exact solution for the asymptotic value of S (as Np -► oo) and 
show that it is a function of the and Np. Close 
examination of Figure 8 reveals the similar form of all three 
curves. That is, for any value of Np the difference between 
^oh,. and (for points) or (Sob,. and 6^,^^ (for particles) 
appears to be the same. Based on the previous observation 
that this difference is due to the particles positioned along 
the border, the following relationship was defined: 

^A6 = 6^^,. - (4.1) 

^ p 

where Npp is the number of particles along the border and A<5 
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AVERAGE FIRST NEAREST NEIGHBOR DISTANCE 



0.10 



0.09 



0.08 



Border particles included in 5 calculations: 
-^Particles, a^ = 0.40, rp = 0.05 
-♦-Particles, a^ = 0 . 10 , rp = 0.05 
—*— Calculated Points 




Border particles excluded from 5 calculations: 



—D— Particles, = 0.40, rp = 0.50 
—0— Particles, a£ = 0.10, rp = 0.05 
—Miscalculated Points 




10 100 1000 10000 

TOTAL NUMBER OF PARTICLES 

Figure 10. Total number of particles vs average first nearest 
neighbor distance for Na “ IvO, The three pairs of curves 
demonstrate the affect of particles located on the border. 
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is the overestimate in S associated with each particle along 
the border. For a random distribution of particles in a 
square region, Np can be taken to be equal to the square of 
the number of particles along the border, Np^^. Therefore Npt, 
= Np'^^ and equation (4.1) can be expressed as: 

— ^A6 = ^obs. ~ ^Theo. (4.2) 

Figure 11 is a plot of Np''^^ vs (Sobs. “ '^Theo. from the data 
compiled in Table 1. These curves clearly define a linear 
relationship from which the following equations can be 
derived: 



Points: = 0.0235^ + (4.3) 

v/^ 

= 0.1: = 0.0248—^ + 0.0100 + (4.4) 

a^ = 0.4: = 0.0248— + 0.0309 + (4.5) 

v/^ 

By taking the limit as Np -► <» in each equation (4.3), (4.4) 
and (4.5), (Sobs. - ^Asym.* The solutions for in these 

equations are 0.0500, 0.0600 and 0.0809, respectively. In 

general, this value for (SAsym. taken to be the average 

first nearest neighbor distance for any distribution of points 
or particles of finite size. 
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- 0-5 

Figure 11. N_ vs 
O.IO and 0.40. 
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5 for area fractions of 
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Based on the results above all subsequent simulations were 
made using both the border particle exclusion algorithm and a 
nominal value of 5000 particles or points. 

C. RANDOM DISTRIBUTION OP CONSTANT SIZE PARTICLES 

For all values of evaluated, 1.275 < < 63.82, the 
calculated distribution for random points very closely 
approximated the corresponding Poisson distribution (Figure 12 
-16). For area fractions less than 0.05, as shown in Figure 
12 , the particle distributions generated approximated the 
Poisson as well. In these dilute distributions the particles 
behave like points since the likelihood of two randomly 
generated particles overlapping is not significant. The 
difference between the distributions is that the particle 
distribution becomes truncated at a value equal to twice the 
particle radius. This was a result of the non-overlap 
criterion preventing the particle centers from getting closer 
than one particle diameter. Figures 17 - 21 are the 
corresponding plots of the actual particle distributions 
represented in Figures 12 - 16. 

Figure 22 demonstrates the effect that, as the area 
fraction is increased, the distribution is shifted towards 
lower values of NND. A shift in S from the theoretical 
Poisson value (^Theo. = towards the value for a 
hexagonal array (6 hcx. = 1.075Na’^^) was also observed. At high 
area fractions the distributions for particles approach delta 
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RELATIVE NUMBER OF PARTICLES 



0.50 



— Poisson Dist., = 1.275 
♦-Calculated Point Dist., = 1.275 
^Const. Size Particles, = 1.273 
a^ = 0.01, Ip = 0.05, 6 = 0.4608 



0.20 — 




0.00 0.50 1.00 1.50 

FIRST NEAREST NEIGHBOR DISTANCE 

Figure 12. Random distribution of points and constant size 
particles. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 



0.30 



0.20 



0.10 



0.00 



Figure 

particl 



— Poisson Dist., N^= 6.367 

Calculated Point Dist., N 6.367 
4 — Const. Size Particles, N^= 6.377 
a^ = 0.05, tp = 0.05, 5 = 0.2214 







0.00 0.50 1.00 

FIRST NEAREST NEIGHBOR DISTANCE 

13. Random distribution of points and constant size 
es. 



1.50 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 



— Poisson Dist., N^= 12.77 
♦-Calculated Point Dist., = 12.77 
4— Const. Size Particles, = 12.72 

a^ = 0.10, r = 0.05, 6 - 0.1699 



0.30 — 




0.00 0.10 0.20 0.30 0.40 0.50 

FIRST NEAREST NEIGHBOR DISTANCE 

Figure 14. Random distribution oi points and constant size 
particles. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 



♦ 



Poisson Dist., N ^ = 38.14 
Calculated Point Dist., N 
Const. Size Particles, N ^ 
a^ = 0.30, Tp = 0.05, 6 ” 



38.14 
= 37.95 
0.1218 



0.30 — 




0.00 0.10 0.20 0.30 0.40 0.50 

FIRST NEAREST NEIGHBOR DISTANCE 

Figure 15. Random distribution of points and constant size 
particles. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 — 



— Poisson Dist., N^= 63.82 
-•-Calculated Point Dist., 

-♦—Const. Size Particles, = 62.82 



= 63.82 




0.00 0.10 0.20 0.30 0.40 0.50 

FIRST NEAREST NEIGHBOR DISTANCE 

Figure 16. Random distribution of points and constant size 
particles. 
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PARTICLE AREA FRACTION: 0.01 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.4608 
NUMBER PARTICLES PER UNIT AREA: 1.275 



Figure 17. Random distribution of constant size particles. 
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PARTICLE AREA FRACTION: 0.05 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.2214 
NUMBER PARTICLES PER UNIT AREA: 6.377 



Figure 18. Random distribution of constant size particles. 
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PARTICLE AREA FRACTION: 0.10 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.1699 
NUMBER PARTICLES PER UNIT AREA: 12.72 



Figure 19. Random distribution of constant size particles. 
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PARTICLE AREA FRACTION: 0.30 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.1218 
NUMBER PARTICLES PER UNIT AREA: 37.95 



Figure 20. Random distribution of constant size particles. 
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PARTICLE AREA FRACTION: 0.50 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.1082 
NUMBER PARTICLES PER UNIT AREA: 62.82 



Figure 21. Random distribution of constant size particles. 
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functions with the minimum first NND having the highest 
frequency of occurrence. 

The effects of increasing area fraction for finite sized 
particles is shown in Figure 23. As increases from values 
of very diluted arrays to more dense distributions the 
corresponding values of <5 shift from that of a random 
distribution of infinitesimal points to that of a hexagonal 
array. Figure 23 shows that this trend is consistent for 
various particle radii. The reader should note the asymptotic 
approach towards the random theoretical distribution of 
points. This is indicative of using a sufficient number of 
particles in the simulations. 

D. NON-RANDOM DISTRIBUTION OF CONSTANT SIZE PARTICLES 

A non-random distribution was simulated by dividing the 
total defined area in to five sub-areas of equal size. In all 
simulations there were three low density zones separated by 
two high density zones. The high density zones were each 
assigned 35% of the particles occupying the total area, while 
the low density zones were each assigned 10%. These 
parameters were used for all simulations. 

The non-random distribution of constant size particles for 
various area fractions are shown in Figures 24 - 27. Plots of 
the particles themselves are shown in Figures 28 - 31. These 
figures were compared with the corresponding distributions and 
particle plots, with respect to area fraction, in the random 
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AVERAGE FIRST NEARST NEIGHBOR DISTANCE 




NUMBER OF PARTICLES PER UNIT AREA 

Figure 23. Number of particles per unit area vs average first 
nearest neighbor distance for particles of varying radius. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 



Poisson Dist., N^= 1.273 

-^Calculated Point Dist., N^ = 1.273 
— 4 — Const. Size Particles, N * 1.273 
a^ = 0.01, tp * 0.05, 5 = 0.4403 

Note: Banded distribution simulated 
by dividing total area into 5 
sub*areas, each containing 10%, 35%, 
10%, 35% and 10% respectively, of 
the total number of points or of the 
total area of particles. 




0.00 0.50 1.00 1.50 

FIRST NEAREST NEIGHBOR DISTANCE 



Figure 24. Non-random distribution of points and constant size 
particles. 
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RELATIVE NUMBER OF PARTICLES 




0.00 0.50 1.00 1.50 

FIRST NEAREST NEIGHBOR DISTANCE 



Figure 25. Non-random distribution of points and constant size 
particles. 
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RELATIVE NUMBER OF PARTICLES 




0.00 0.10 0.20 0.30 0.40 0.50 

FIRST NEAREST NEIGHBOR DISTANCE 



Figure 26. Non>random distribution of points and constant size 
particles. 
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RELATIVE NUMBER OF PARTICLES 



0.70 




FIRST NEAREST NEIGHBOR DISTANCE 



Figure 27. Non-random distribution oi points and constant size 
particles. 
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TOTAL PARTICLE AREA FRACTION: 0.01 
SUB-AREA PARTICLE AREA FRACTIONS: 0.10, 0.35, 
0.10, 0.35, 0.10 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.4403 
NUMBER PARTICLES PER UNIT AREA: 1.273 



Figure 28. Non— random distribution of constant size particles. 
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TOTAL PARTICLE AREA FRACTION: 0.05 
SUB-AREA PARTICLE AREA FRACTIONS: 0.10, 0.35, 
0.10, 0.35, 0.10 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.2110 
NUMBER PARTICLES PER UNIT AREA: 6.383 



Figure 29. Non-random distribution of constant size particles. 
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TOTAL PARTICLE AREA FRACTION: 0.10 
SUB-AREA PARTICLE AREA FRACTIONS: 0.10, 0.35, 
0.10, 0.35, 0.10 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.1640 
NUMBER PARTICLES PER UNIT AREA: 12.78 



Figure 30. Non— random distribution of constant size particles. 
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TOTAL PARTICLE AREA FRACTION: 0.30 
SUB-AREA PARTICLE AREA FRACTIONS: 0.10, 0.35, 
0.10, 0.35, 0.10 
PARTICLE RADIUS: 0.05 

AVERAGE NEAREST NEIGHBOR SPACING: 0.1187 
NUMBER PARTICLES PER UNIT AREA: 38.14 



Figure 31. Non-random distribution of constant size particles. 
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distribution case. It was observed that the fundamental form 
of the distributions, trend in values of the observed S, and 
first NND were similar. 

Additionally, in the non-random case, there was a small 
but consistent shift in the distribution towards lower values 
of first NND. This pattern was attributed to the presence of 
the high density bands. For a given N^, the greater degree of 
packing of points or particles in the high density zones 
caused the first NND to shift toward lower values. This 
effect was also reflected in a greater frequency of the 
minimum first NND when compared to that of the random case. 
A summary of the observed trends in the non-random/ constant 
particle size distributions are shown in Figure 32. 

The corresponding particle plots, comparing the relative 
densities of the high and low density bands, are shown in 
Figures 28 - 31. Compare these plots with those of the random 
distributions for the same area fractions (Figures 17 - 21) . 
For af = 0.3 0 the difference between the random (Figure 20) and 
non-random (Figure 31) distributions is completely obvious to 
the eye. There is a clear distinction between the high and 
low density bands in the non-random distribution. This 
distinction makes it easy for the observer to tell the 
difference between the random and non-random particle 
distributions based on these plots. Again, for a^ = 0.30, 
compare the distribution curves for the random (Figure. 15) 
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FIRST NEAREST NEIGHBOR DISTANCE 



and non-random (Figure 27) cases. The distinction is not 
great and the discernable difference is a greater maximum 
frequency for the lowest first NND in the non-random case. 

The visual difference in the particle pldts is much more 
subtle for low area fractions compared to high area fractions. 
Compare the random and non-random particle plots for af = 0.01 
(Figures 17 and 28) . Unlike the plots for the 30% area 
fraction, the distinction is not obvious. However, the same 
can be said for the NND distribution curves (Figures 12 and 
24). For the low area fraction both curves approach the 
distribution of random points. This made any distinction 
between them even more subtle than the curves for the 30% area 
fraction. 

E. RANDOM DISTRIBUTION OF PARTICLES WITH A LOGNORMAL SIZE 

DISTRIBUTION 

In actual microstructures the second phase-particle size 
generally will not be constant throughout the material. In 
fact it has often been reported that second-phase particles 
follow a lognormal size distribution. Such was the basis for 
the following experiments. 

The same area fractions were analyzed for particles with 
a lognormal size distribution (Figures 33 - 37) as with those 
for constant size particles. The corresponding particle plots 
are shown in Figures 38 - 42. For very dilute arrays of 
particles, the distributions approached very closely to those 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 



— Poisson Disl., N^ = 1.275 
•—Calculated Point Dist., N^= 1.275 
-♦—Lognormal Size Dist., = 1.273 
a, = 0.01, r = 0.0499, 6 = 0.2214 

o - 0.4 



0.30 



0,20 — 



0.10 



0.00 




0.00 0.50 1.00 1.50 

FIRST NEAREST NEIGHBOR DISTANCE 



Figure 33. Random distribution of points and particles with a 
lognormal size distribution. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 



— Poisson Dist., N^ = 6.367 
•—Calculated Point Dist., N^= 6.367 
^Lognormal Size Dist., = 5-356 
a^ = 0.05, = 0.0505, 6 = 0.2382 

a - 0.4 



0.30 — 




0.00 0.50 1.00 1.50 

FIRST NEAREST NEIGHBOR DISTANCE 



Figure 34. Random distribution of points and particles with a 
lognormal size distribution. 



58 



RELATIVE NUMBER OF PARTICLES 




0.00 0.10 0.20 0.30 0.40 0.50 

FIRST NEAREST NEIGHBOR DISTANCE 



Figure 35. Random distribution of points and particles with a 
lognormal size distribution. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 



— Poisson Dist., N^ = 38.14 
*- Calculated Point Dist., N^= 38.14 
Lognormal Size Dist., = 31.56 
a^ = 0.30, = 0.0506, 6 - 0.1256 

a - 0.4 



0.30 




0.00 0.10 0.20 0.30 0.40 0.50 

FIRST NEAREST NEIGHBOR DISTANCE 



Figure 36. Random distribution of points and particles with a 
lognormal size distribution. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



Poisson Dist., N^ = 63.82 

—■—Calculated Point Dist., N^= 63.82 
—♦—Lognormal Size Dist., = 59.65 




0.00 0.10 0.20 0.30 0.40 0.50 



FIRST NEAREST NEIGHBOR DISTANCE 



Figure 37. Random distribution oi points and particles with a 
lognormal size distribution. 
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PARTICLE AREA FRACTION: 0.01 
MEAN PARTICLE RADIUS: 0.0499 
AVERAGE NEAREST NEIGHBOR SPACING: 0.4944 
NUMBER PARTICLES PER UNIT AREA: 1.094 

Figure 38. Random distribution of particles with a lognormal 
size distribution. 
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PARTICLE AREA FRACTION: 0.05 
MEAN PARTICLE RADIUS: 0.0505 
AVERAGE NEAREST NEIGHBOR SPACING: 0.2382 
NUMBER PARTICLES PER UNIT AREA: 5.356 

Figure 39. Random distribution of particles with a lognormal 
size distribution. 
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PARTICLE AREA FRACTION: 0.10 
MEAN PARTICLE RADIUS: 0.0503 
AVERAGE NEAREST NEIGHBOR SPACING: 0.1808 
NUMBER PARTICLES PER UNIT AREA: 10.75 

Figure 40. Random distribution of particles with a lognormal 
size distribution. 
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PARTICLE AREA FRACTION: 0.30 
MEAN PARTICLE RADIUS: 0.0506 
AVERAGE NEAREST NEIGHBOR SPACING; 0.1256 
NUMBER PARTICLES PER UNIT AREA: 31.56 

Figure 41. Random distribution of particles with a lognormal 
size distribution. 
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PARTICLE AREA FRACTION: 0.50 
MEAN PARTICLE RADIUS: 0.0479 
AVERAGE NEAREST NEIGHBOR SPACING: 0.1012 
NUMBER PARTICLES PER UNIT AREA: 59.65 

Figure 42. Random distribution of particles with a lognormal 
size distribution. 
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calculated for points and to the Poisson distribution (Figure 
33) . But now, as the area fraction is increased the particle 
distribution still maintains a shape similar to the Poisson 
distribution. Unlike the case for constant size particles, 
the minimum first NND can be smaller than the mean particle 
diameter. A smaller particle can always be fitted between two 
or more larger ones, thus reducing the first NND below the 
value of the mean particle diameter. This also resulted in a 
more nearly constant peak frequency of first NND for all 
ranges of area fractions examined. 

The trends for 6 of the particles and ranges of first NND 
(Figure 43) were the same as those in the constant size 
particle case for lower values of N^. At larger values of N^ 
the value of S for lognormally distributed particles becomes 
less than that for the corresponding constant particle sizes. 

F. NON-RANDOM DISTRIBUTION OF PARTICLES WITH A LOGNORMAL SIZE 

DISTRIBUTION 

The same parameters for banding used in the non-random, 
constant size particle case were used in the case of non- 
random distribution of particles with a lognormal size 
distribution. The results of this analysis are shown in 
Figures 44 - 47 for the distributions and Figures 48 - 51 for 
the particle plots. The features are similar in many respects 
to those reported earlier for non-random/constant particle 
size distributions. As the area fraction increases S shifts 
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FIRST NEAREST NEIGHBOR DISTANCE 



RELATIVE NUMBER OF PARTICLES 




FIRST NEAREST NEIGHBOR DISTANCE 



Figure 44. Non*random distribution of points and particles with a 
lognormal size distribution. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



0.40 



Poisson Dist., N^ = 6.378 

—•—Calculated Point Dist., = 6.378 
—♦—Lognormal Size Dist., = 5.377 

^’m " 0.0503, 6 - 0.2290 

a - 0.4 

Note: Banded distribution simulated by 
dividing the total area into 5 sub-areas, 
each containing 10%, 35%, 10%, 35% and 
10% respectively, of the total number of 
points or of the total area of particles. 
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0.50 



1.00 
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Figure 45. Non-random distribution of points and particles 
lognormal size distribution. 
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RELATIVE NUMBER OF PARTICLES 



0.50 



Poisson Dist., Na = 12.77 
Calculated Point Dist., = 12.77 



-♦—Lognormal Size Dist., 



= 10.65 



0.40 — 



a^ = 0.10, = 0.0505, 

5 - 0 . 1745 , o - 0.4 

Note: Banded distribution 
simulated by dividing the total 
area into 5 sub-areas, each 
containing 10%, 35%, 10%, 35% 
and 10% respectively, of the 
total number of points or of 
the total area of particles. 
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Figure 46. Non-random distribution of points and particles with a 
lognormal size distribution. 
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Figure 
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— Poisson Dist., N^= 38.31 
•—Calculated Point Dist., = 38.31 
-♦—Lognormal Size Dist., = 37.79 
a^ = 0.30, = 0.0470, 6 -O.llOO 

a = 0.4 
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47. Non-random distribution of points and particles 
mal size distribution. 
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TOTAL PARTICLE AREA FRACTION: 0.01 
SUB-AREA PARTICLE AREA FRACTIONS; 0.10, 0.35, 

0.10, 0.35, 0.10 

MEAN PARTICLE RADIUS: 0.0502 

AVERAGE NEAREST NEIGHBOR SPACING: 0.4703 

NUMBER PARTICLES PER UNIT AREA: 1.074 

Figure 48. Non— random distribution of particles with a log- 
normal size distribution. 
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TOTAL PARTICLE AREA FRACTION: 0.05 
SUB-AREA PARTICLE AREA FRACTIONS: 0.10, 0.35, 

0.10, 0.35, 0.10 

MEAN PARTICLE RADIUS: 0.0503 

AVERAGE NEAREST NEIGHBOR SPACING: 0.2290 

NUMBER PARTICLES PER UNIT AREA: 5.377 

Figure 49. Non— random dislribuiion of particles with a log- 
normal size distribution. 



74 




TOTAL PARTICLE AREA FRACTION: 0.10 
SUB-AREA PARTICLE AREA FRACTIONS: 0.10, 0.35, 

0.10, 0.35, 0.10 

MEAN PARTICLE RADIUS: 0.0508 

AVERAGE NEAREST NEIGHBOR SPACING: 0.1762 

NUMBER PARTICLES PER UNIT AREA: 10.53 

Figure 50. Non-random distribution of particles with a log- 
normal size distribution. 
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TOTAL PARTICLE AREA FRACTION: 0.30 
SUB-AREA PARTICLE AREA FRACTIONS: 0.10, 0.35, 
0.10, 0.35, 0.10 

MEAN PARTICLE RADIUS: 0.0470 

AVERAGE NEAREST NEIGHBOR SPACING: 0.1100 

NUMBER PARTICLES PER UNIT AREA: 37.79 

Figure 51. Non— random distribution of particles with a log- 
normal size distribution. 
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from a very close approximation to the value for a Poisson 
distribution towards that of a hexagonal array. Unlike the 
cases for constant size particles, the distribution curve 
maintains a shape similar to the Poisson distribution and does 
not reach a peak frequency of first NND corresponding to the 
minimum first NND. This, again, is due to the possibility of 
positioning a smaller particle between two or more larger 
ones. There is a consistent increase in the maximum 
occurrence of any first NND value over that of the random 
case. This is also a result of the greater density of 
particles in the high density zones. Figure 52 summarizes the 
shift in the distributions as the area fraction increases. 

In comparing the particle plots for non-random 
distributions of particles with a lognormal size distribution. 



Figures 48 - 51, 


to 


those of the random 


case, the 


same 


observations are 


made 


as for the constant 


particle 


size 


distributions . 


For 


relatively high area 


fractions 


the 



difference between random and non-random particle dispersions 
is obvious to the eye. But for very low area fractions it 
becomes more difficult to differentiate between the two. This 
observation is summarized graphically in Figures 53, 54 and 
55. The first two figures represent random and non-random 
distributions, respectively, for particles of constant and 
varying size. Both figures are indeed very similar, each 
showing the effect of smaller particles positioned between 
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0.00 0.50 1.00 1.50 2.00 2.50 

FIRST NEAREST NEIGHBOR DISTANCE 

Figure 52. Non-random particle distributions for various area fractions; 
lognormal particle size distribution; mean particle radius “ 0.05; standard 
deviation ■ 0.4. 



3 -a AVERAGE FIRST NEAREST NEIGHBOR DISTANCE 



1 



•—Particle Radius = 0.05 
-c^Lognormal Particle Size Dist., 
rjj^ = 0.05, a = 0.40 




1 10 100 

NUMBER OF PARTICLES PER UNIT AREA 

igure 53. Number of particles per unit area vs average first 
earest neighbor distance; for random distributions of constant 
size particles and for particles with a lognormal size distribution. 
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AVERAGE FIRST NEAREST NEIGHBOR DISTANCE 




NUMBER OF PARTICLES PER UNIT AREA 

Figure 54. Number of particles per unit area vs average first 
nearest neighbor distance; for non-random distributions of constant 
size particles and for particles with a lognormal size distribution. 
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AVERAGE FIRST NEAREST NEIGHBOR DISTANCE 



1 



« — Random Distribution 
-* — Non-random Distribution 




1 10 100 

NUMBER OF PARTICLES PER UNIT AREA 

Figure 55. Number of particles per unit area vs average first 
nearest neighbor distance; for random and non-random 
distributions of particles with a log normal size distribution, 
tm “ 0-05, Q - 0.4. 
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larger ones reducing the value of S for a given N^. Further, 
in comparing the two figures the minimal deviation between 
random and non-random distributions was reflected in the close 
similarities of the curves. This observation is more apparent 
in figure 55 for the case of random and non-random 
distributions of particles with a lognormal size variance. 

The analysis conducted shows that the human eye is very 
capable of detecting non-randomness of particle distributions 
for sufficiently dense arrays. However, when presented with 
the statistical data, the discernment of random or non-random 
is not so obvious . The degree of non-randomness imposed on 
the model was judged to be relatively high, 35% and 10% for 
high and low density bands respectively. Had the difference 
between the two bands been less, there would have been even 
greater similarities between the random and non-random 
distributions. The question then becomes; is it possible to 
develop a system which can differentiate between random and 
non-random particle distributions and quantify the results? 
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V. CONCLUSIONS AND RECOMMENDATIONS 

The following conclusions may be drawn from the work 
described in the preceding sections. 

1. The simulations performed for large values of 
reproduce the work performed by Pas [Ref. 5] and therefore 
substantiate those results. 

2. For more dilute arrays, because the problem of 
particles along the border was not recognized, Pas's 
calculations for first NND was overestimated. 

3. As the number of particles utilized in any simulation 

increases, the relative number of particles appearing along 
the border decreases. For this reason the errors resulting in 
calculations of first NND can be reduced by using a sufficient 
number of particles (at least 1000) . Calculations of first 
NND for points approached the theoretical Poisson distribution 
as the total number of points utilized was increased. The 
same can be said for the distribution of finite sized 
particles. However in this case, as the total number of 

particles increases (approaches 1000 and beyond) , S 
asymptotically approaches a constant value which lies between 
the minimum possible spacing (2 x tp) and the spacing for a 
hexagonal array. Based on these observations it is possible 
to compensate for the statistical inaccuracies by utilizing at 
least 1000 particles in the simulations. 
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4. Attempts to compensate for errors arising in the first 
NND calculations when using a small number of particles by 
eliminating particles along the border was only partially 
successful. Any gains made using the procedure to reduce the 
border particle problem for small number of particles was 
overtaken by simply increasing the total number of particles. 

5. In all four cases analyzed (random 
distribution/constant size particles, non-random 
d i s t r i but i on / con s t a nt size particles, random 
distribution/ lognormal particle size distribution, and non- 
random distribution/ lognormal particle size distribution) the 
non-overlap criterion caused the distribution to deviate from 
the Poisson to that of a hexagonal array as the density of the 
particles increased. Using a nominal value of 5000 particles 
the model very closely simulated the Poisson distribution in 
all cases for dilute arrays. 

6. For constant size particles in high area fractions (af 
> 0.10), for both random and non-random distributions, the 
most frequently occurring first NND was the minimum first NND. 
This was due to the non-overlap criterion and that the minimum 
first NND can not be less than the particle diameter. 

7. The non-random distribution of points and constant 
size particles was shifted to the left of the corresponding 
Poisson distribution. This was due to the presence of the 
high density zones. For a given N,^, the closer packing of the 
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points or particles in the high density zones causes NND to 
Ishift towards lower values. This was also reflected in a 
igreater frequency of the minimum first NND when compared to 
that of the random case. 

8. The random distribution of particles with a lognormal 

i 

size distribution maintained the shape of a Poisson field for 
all area fractions. Unlike the case for constant size 
particles, the minimum first NND can be smaller than the mean 
particle diameter. A smaller particle can always fit between 
two or more larger ones thus lowering the minimum first NND. 
This also resulted in a more constant peak frequency of first 
NND for all ranges of area fractions examined. These 
observations held true as well for the case of non- 
random/lognormal particle size distribution. 

After reviewing all the accumulated data there did not 
appear to be as significant a difference between the plotted 
distributions of the random and non-random cases as there was 

I 

between the plots of the particles themselves. Based then on 
the statistical nature of this work it appears that the best 
dif ferentiater between random and non-random particle 
distributions is the human eye. 

To follow up on this work the following studies are 
recommended. A similar analysis should be made for the three 
dimensional case in order to parallel the actual particle 
distributions. Other types of non-randomness, such as 
clusters, should be explored to determine their effects on 



85 



distributions. Area fractions greater than 0.50 should be 
examined as well. Finally, an analysis should be made, based 
on results of this work, to determine the feasibility of 
developing a computer program which could determine the degree 
of randomness in the distributions of second phase particles. 
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APPENDIX A. PROGRAM RANPARTl.FOR 





PROGRAM 


RANPARTl 


* 


Glossary 


of variables: 


* 


Pi 


3.14159 


* 


Xpos 


X coordinate of an interior particle 


* 


Ypos 


Y coordinate of an interior particle 


* 


XposB 


X coordinate of a particle along the 


* 




border 


* 


YposB 


Y coordinate of a particle along the 


* 




border 


ic 


AP 


Area percentage of the particles 


* 


TPA 


Total area fraction of the particles 


* 


Pspace 


Min. spacing between particles 


* 


Pradius 


Particle radius 


* 


TOTarea 


Defined area where particles are positioned 


* 


Parea 


Area occupied by each particle 


* 


dp 


Distance between interior particles 


* 


dp 2 


Distance between interior particles 


* 


dp2B 


Distance between an interior particle and a 


* 




border particle 


* 


dp 3 


Nearest neighbor spacing 


* 


dp 4 


Sum of nearest neighbor spacings 


* 


dist 


Particle diameter plus min spacing 


* 


dpavg 


Average nearest neighbor spacing 


* 


X 


Horizontal dimension of defined area 


* 


Y 


Vertical dimension of defined area 


* 


dpmin 


Minimum nearest neighbor spacing 


* 


dp2Bmin 


Distance between an interior particle and 


* 




its nearest border particle 


* 


dpmax 


Maximum nearest neighbor spacing 


* 


dpinc 


Size of each range between dpmin and dpmax 


* 




based on the Sturgis Rule 


* 


Delta 


Minimum value of each range in dpinc 


* 


Np 


Number of particles corresponding to each 


* 




Delta 


* 


NPA 


Number of particles per unit area 


* 


Pnum 


Total number of particles generated 


* 


Ans 


Either "A" for area or "p" for points. 


* 




depending on 


* 




type of particles 


* 


Part 


Either ”n" for a preset number of particles 


* 




or "r” for a preset particle radius 



* 



REAL Pi, Xpos(5200), Ypos(5200), XposB(250), YposB(250), 
:TPA, AP, Pspace, Pradius, TOTarea, Parea, dp, dp2 , dp2B, 
:dp2Bmin, dp3(5200), dp4, dist, dpavg, X, Y, dpmin, NPA, 
:dpmax, dpinc, g, Delta(30) , Np(30) 
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INTEGER a, b, c, d, i, j, k, 1, m, n, Pnum, f, h, e 
CHARACTER Ans*l, Part*l 
PARAMETER (Pi=3. 14159) 

DATA n, dp4, e, dpmin, dpmax, Np(l) , a, b / 

: 1 , 0 . 0 , 0 , 10 . 0 , 0 . 001 , 0 , 1 , 1 / 

OPEN (10, f ile='rplx.dat' , status='new' ) 

OPEN (20, f ile=' rply.dat ' , status='new' ) 

OPEN (30 , f ile='Dell. dat ' , status='new' ) 

OPEN (40, f ile='Npl .dat ' , status='new' ) 

OPEN (50, f ile='rpl .dat' , status='new' ) 

PRINT*,' Enter the area percentage of particles.' 
READ*, AP 

PRINT*,' Enter the length of the X - scale.' 

READ*, X 

PRINT*,' Enter the length of the Y - scale.' 

READ*, Y 

TOTarea = X*Y 

TPA = AP/100*TOTarea 

PRINT*,' Enter the minimum spacing between particles.' 
READ*, Pspace 

PRINT*,' If particle generation is to based on a preset 
: number of' 

PRINT*,' particles enter "n” . OR If particle 
: generation is to' 

PRINT*,' be based on a preset particle radius enter 
READ*, Part 

IF (Part .EQ. 'n') THEN 

PRINT* , 'Enter total number of particles.' 

READ*, Pnum 
Parea = TPA/Pnum 
Pradius = SQRT(Parea/Pi) 

ELSEIF (Part .EQ. 'r') THEN 

PRINT*,' Enter the radius of the particles.' 

READ*, Pradius 

Parea = Pi* (Pradius**2) 

Pnum = NINT (TPA/Parea) 

END IF 

PRINT*, ' Are the particles finite areas or points? 

: Enter "A” or ”p” . ' 

READ*, Ans 

Determine random position of particles. 

PRINT*,' Determining random position of particles.' 
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dist = 2*Pradius + Pspace 
Xpos(l) = RRANDO 
DO 100 i = 1, Pnum 
150 CONTINUE 

Xpos(i) = RND()*X 
Ypos(i) = RND() *Y 

* 

* Check for overlap of particles with a finite area 

* 

IF (Ans .EQ. 'p') GOTO 250 
DO 200 k = 1, i - 1 

dp = SQRT( (Xpos(i) - Xpos(k))**2 + 

: (Ypos(i) - Ypos(k))**2) 

IF (dp .LT. dist) GOTO 150 
200 CONTINUE 

250 CONTINUE 

100 CONTINUE 
* 

* Segregate particles along the perimeter from those in 

* the interior of the specified area. 

* 

PRINT*,' Segregating particles on the border from the 
: interior. ' 

DO 275 d = 1, Pnum 

IF ((Xpos(d) .LE. Pradius + Pspace) .OR. 

: (Xpos(d) .GE. X - Pradius - Pspace) .OR. 

: (Ypos(d) .LE. Pradius + Pspace) .OR. 

: (Ypos(d) .GE. Y - Pradius - Pspace)) THEN 

XposB(b) = Xpos(d) 

YposB(b) = Ypos(d) 
b = b + 1 
ELSE 

Xpos(a) = Xpos(d) 

Ypos(a) = Ypos(d) 

WRITE (10,*) Xpos(a) 

WRITE(20,*) Ypos(a) 
a = a + 1 
ENDIF 

275 CONTINUE 
* 

* Determine the average distance between nearest 

* neighbors. 

* 

PRINT*,' Determining nearest neighbor distances.' 

DO 300 j = 1, Pnum - (b - 1) 
dp3(n) = SQRT(X**2 + Y**2) 
dp2Bmin = SQRT(X**2 + Y**2) 

DO 350 c = 1, (b - 1) 

dp2B = SQRT( (Xpos(j) - XposB(c))**2 + 

: (Ypos(j) - YposB(c) ) **2) 

dp2Bmin = MIN(dp2Bmin, dp2B) 
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350 CONTINUE 

DO 400 1 - 1 , Pnum - ( b - 1 ) 

dp2 = SQRT( (Xpos(l) - Xpos(j))**2 + 

: (Ypos(l) - Ypos(j))**2) 

IF (dp2 .NE. 0.0) THEN 

dp3(n) = MIN(dp2, dp3(n), dp2Binin) 

ENDIF 

400 CONTINUE 

n = n + 1 
300 CONTINUE 

DO 500 m = 1, Pnum - (b - 1) 
dp4 = dp4 + dp3 (m) 

500 CONTINUE 

dpavg = dp4/(Pnum - (b - 1)) 

* 

* Determine the number of particles falling in each range 

* of nearest neighbor spacings (based on the Sturgis 

* rule) . 

* 

PRINT*,' Grouping nearest neighbor values into ranges.' 
DO 600 f = 1, Pnum - 1 - (b - 1) 

dpmin = MIN(dp3(f), dp3(f+l), dpmin) 
dpmax = MAX(dp3(f), dp3(f+l), dpmax) 

600 CONTINUE 

dpinc = (dpmax - dpmin) /(I + 3 . 3*LOG10 (Pnum - (b - 1) ) ) 
DO 700 g = dpmin, dpmax, dpinc 
e = e + 1 

DO 800 h = 1, Pnum - (b - 1) 

IF ((dp3(h) .GE. g) .AND. (dp3 (h) .LT. g + dpinc)) 

: THEN 

Np(e) = Np(e) + 1 
ENDIF 

800 CONTINUE 

Delta(e) = (2*g + dpinc) /2 
WRITE ( 30,*) Delta (e) 

WRITE(40,*)Np(e)/(Pnum - (b - 1) ) 

700 CONTINUE 

NPA = (Pnum - (b - 1))/(X - 2*Pradius - 2*Pspace)/ 

: (Y - 2*Pradius - 2*Pspace) 

WRITE (50, 10) 'The area percentage of particles is', AP 
WRITE(50,20) 'X dimension is',X 
WRITE(50,20) 'Y dimension is',Y 

WRITE(50,20) 'Min spacing between particles is',Pspace 
WRITE(50, 20) 'Particle radius is', Pradius 
WRITE (50, 20) 'Particle area is', Parea 
WRITE(50, 20) 'Total number of particles is', Pnum 
WRITE (50, 20) 'Number of particles along the perimeter 
: is ' , b - 1 

WRITE (50, 20) 'Number of particles in area interior is', a 
1 

WRITE (50 , 20) 'Number particles per unit area', NPA 
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10 

20 



WRITE (50, 20) 'Average spacing between nearest neighbors 
:is', dpavg 



WRITE(50, 20) 'Min nearest neighbor spacing is',dpmin 
WRITE(50, 20) 'Max nearest neighbor spacing is',dpmax 
PRINT*,' The area percentage of particles is', AP 

PRINT*,' X dimension is',X 

PRINT*,' Y dimension is',Y 

PRINT*,' Min spacing between particles is',Pspace 

PRINT*,' Particle radius is', Pradius 
PRINT*,' Particle area is', Parea 
PRINT*,' Total number of particles is', Pnum 
PRINT*,' Number of particles along the perimeter', b - 1 
PRINT*,' Number of particles in area interior', a - 1 
PRINT*,' Number particles per unit area is', NPA 
PRINT*,' Average spacing between nearest neighbors is'. 



: dpavg 

PRINT*,' Min nearest neighbor spacing is',dpmin 
PRINT*,' Max nearest neighbor spacing is',dpmax 
FORMAT ( //////,T8,A,F10.4) 

FORMAT (T8,A,F10.4) 

END 
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APPENDIX B 



PROGRAM RANPART2.FOR 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



PROGRAM RANPART2 



Glossary of variables: 

Pi 3.14159 

Xpos X coordinate of an interior particle 

Ypos Y coordinate of an interior particle 

XposB X coordinate of a particle along the border 

YposB Y coordinate of a particle along the border 

AP Area percentage of the particles 

TPA Total area fraction of the particles 

Pspace Min. spacing between particles 

Pradius Particle radius 

TOTarea Defined area of particles 

Parea Area occupied by each particle 

dp Distance between interior particles 

dp2 Distance between interior particles 

dp2B Distance between an interior particle and a 

border particle 

dp3 Nearest neighbor spacing 

dp4 Sum of nearest neighbor spacings 

dist Particle diameter plus min spacing 

dpavg Average nearest neighbor spacing 

X Horizontal dimension of total area 

Y Vertical dimension of total area 

dpmin Minimum nearest neighbor spacing 

dp2Bmin Distance between an interior particle and 

its nearest border particle 
dpmax Maximum nearest neighbor spacing 

dpinc Size of each range between dpmin and dpmax 

based on the Sturgis Rule 
Delta Minimum value of each range in dpinc 

Np Number of particles corresponding to each 

Delta 

NPA Number of particles per unit area 

Pnum Total number of particles generated 

Ans Either "A” for area or "p” for points, 

depending on type of particles 
Part Either "n" for a preset number of particles 

or "r" for a preset particle radius 
AFl-5 Area percentage of particles in each of the 

five sub-areas 

Al-5 Subtotal of particles in each sub-area 



REAL Pi, Xpos(5200), Ypos(5200), XposB(250), YposB(250), 
:TPA, AP, Pspace, Pradius, TOTarea, Parea, dp, dp2, dp2B, 
;dp2Bmin, dp3(5200), dp4, dist, dpavg, X, Y, dpmin, NPA, 
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:dpmax, dpinc, g, Delta(30), Np(30), Al, AFl, A2 , AF2 , 
:A3, AF3, A4, AF4 , A5, AF5 

INTEGER a, b, c, d, i, j, k, 1, m, n, Pnum, f, h, e 

CHARACTER Ans*l, Part*l 

PARAMETER (Pi=3. 14159) 

DATA n, dp4, e, dpmin, dpmax, Np(l), a, b, Al, A2, A3, 
:A4, A5 / 

: 1 , 0 . 0 , 0 , 10 . 0 , 0 . 001 , 0 , 1 , 1 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 

: 0 . 0 / 

OPEN (10, file='rp2x.dat' , status='new' ) 

OPEN (20, file=' rp2y.dat' , status='new' ) 

OPEN (30, file='Del2 .dat' , status=' new' ) 

OPEN (40,file='Np2.dat' ,status='new') 

OPEN (50, file='rp2 .dat' , status='new' ) 

PRINT*,' Enter the area percentage of particles.' 
READ*, AP 

PRINT*,' Enter 5 sub-area percentages of particles.' 

READ*, AFl, AF2, AF3 , AF4 , AF5 

AFl = AFl/100 

AF2 = AF2/100 

AF3 = AF3/100 

AF4 = AF4/100 

AF5 = AF5/100 

PRINT*,' Enter the length of the X - scale.' 

READ*, X 

PRINT*,' Enter the length of the Y - scale.' 

READ*, Y 

TOTarea = X*Y 

TPA = AP/100*TOTarea 

PRINT*,' Enter the minimum spacing between particles.' 
READ* , Pspace 

PRINT*,' If particle generation is to based on a preset 
: number of' 

PRINT*,' particles enter ”n". OR If particle 
; generation is to' 

PRINT*,' be based on a preset particle radius enter 

• lly*!! / 

• X. • 

READ*, Part 

IF (Part .EQ. 'n') THEN 

PRINT* , 'Enter total number of particles.' 

READ* , Pnum 
Parea = TPA/ Pnum 
Pradius = SQRT( Parea/ Pi) 

ELSEIF (Part .EQ. 'r') THEN 

PRINT*,' Enter the radius of the particles.' 
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READ*, Pradius 

Parea = Pi* (Pradius**2) 

Pnum = HINT (TPA/ Parea) 

ENDIF 

PRINT*,' Are the particles finite areas or points? 

: Enter "A” or "p” . ' 

READ*, Ans 

* 

* Determine random position of particles. 

* 

PRINT*,' Determining random-banded position of 
rparticles. ' 

dist = 2*Pradius + Pspace 
Xpos(l) = RRANDO 
DO 100 i = 1, Pnum 
150 CONTINUE 

Xpos(i) = RND()*X 
Ypos(i) = RND()*Y 

* 

* Check for overlap of particles with a finite area 

* 

IF (Ans .EQ. 'p') GOTO 250 
DO 200 k = 1, i - 1 

dp = SQRT( (Xpos(i) - Xpos(k))**2 + 

: (Ypos(i) - Ypos(k))**2) 

IF (dp .LT. dist) GOTO 150 

200 CONTINUE 

250 CONTINUE 



IF ((Ypos(i) .GE. 0.0) .AND. 


(Ypos(i) .LT. 


0.2*Y) ) 


THEN 


IF (A1 .GT. TPA*AF1) GOTO 
A1 = A1 + Parea 


150 








ELSEIF ((Ypos(i) .GE. 0.2*Y) 
:THEN 


.AND. 


(Ypos(i) 


.LT. 0. 


4*Y)) 


IF (A2 .GT. TPA*AF2) GOTO 
A2 = A2 + Parea 


150 








ELSEIF ((Ypos(i) .GE. 0.4*Y) 
:THEN 


.AND. 


(Ypos(i) 


.LT. 0. 


6*Y) ) 


IF (A3 .GT. TPA*AF3) GOTO 
A3 = A3 + Parea 


150 








ELSEIF ((Ypos(i) .GE. 0.6*Y) 
:THEN 


.AND. 


(Ypos(i) 


.LT. 0. 


8*Y) ) 


IF (A4 .GT. TPA*AF4) GOTO 


150 









A4 = A4 + Parea 

ELSEIF ((Ypos(i) .GE. 0.8) .AND. (Ypos(i) .LE. Y) ) THEN 
IF (A5 .GT. TPA*AF5) GOTO 150 
A5 = A5 + Parea 
ENDIF 

100 CONTINUE 
* 

* Segregate particles along the perimeter from those in 

* the interior of the specified area. 
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PRINT*,' Segregating particles on the border from the 
: interior. ' 

DO 275 d = 1, Pnum 

IF ((Xpos(d) .LE. Pradius + Pspace) .OR. 

: (Xpos(d) .GE. X - Pradius - Pspace) .OR. 

: (Ypos(d) .LE. Pradius + Pspace) .OR. 

: (Ypos(d) .GE. Y - Pradius - Pspace)) THEN 

XposB(b) = Xpos(d) 

YposB(b) = Ypos(d) 

b = b + 1 
ELSE 

Xpos(a) = Xpos(d) 

Ypos(a) = Ypos(d) 

WRITE(10,*) Xpos(a) 

WRITE(20,*) Ypos(a) 
a = a + 1 
ENDIF 

275 CONTINUE 

* 

* Determine the average distance between nearest 

* neighbors. 

* 

PRINT*, ' Determining nearest neighbor distances. ' 

DO 300 j = 1, Pnum - (b - 1) 
dp3(n) = SQRT(X**2 + Y**2) 
dp2Bmin = SQRT(X**2 + Y**2) 

DO 350 c = 1, (b - 1) 

dp2B = SQRT( (Xpos( j) - XposB(c))**2 + 

: (Ypos(j) - YposB(c) ) **2) 

dp2Bmin = MIN(dp2Bmin, dp2B) 

350 CONTINUE 

DO 400 1=1, Pnum - (b - 1) 

dp2 = SQRT( (Xpos(l) - Xpos(j))**2 + 

; (Ypos(l) - Ypos(j))**2) 

IF (dp2 .NE. 0.0) THEN 
dp3(n) = MIN(dp2, dp3(n), dp2Bmin) 

ENDIF 

400 CONTINUE 

n = n + 1 

300 CONTINUE 

DO 500 m = 1, Pnum - (b - 1) 
dp4 = dp4 + dp3 (m) 

500 CONTINUE 

dpavg = dp4/(Pnum - (b - 1)) 

* 

* Determine the number of particles falling in each range 

* of nearest neighbor spacings (based on the Sturgis 

* rule) . 

k 

DO 600 f = 1, Pnum - 1 - (b - 1) 
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600 



800 



700 



dpmin = MIN(dp3(f), dp3(f+l), dpitiin) 
dpmax = MAX(dp3(f), dp3(f+l), dpmax) 

CONTINUE 

dpinc = (dpmax - dpmin) /(I + 3 . 3*LOG10 (Pnum - (b - 1))) 
DO 700 g = dpmin, dpmax, dpinc 
e = e + 1 

DO 800 h = 1, Pnum - (b - 1) 

IF ((dp3(h) .GE. g) .AND. (dp3 (h) .LT. g + dpinc)) 
THEN 



Np(e) = Np(e) + 1 
ENDIF 
CONTINUE 

Delta(e) = (2*g + dpinc) /2 
WRITE (30,*) Delta ( e) 

WRITE (40, *) Np (e) / (Pnum - (b - 1) ) 

CONTINUE 

NPA = (Pnum - (b - 1))/(X - 2*Pradius - 2*Pspace)/ 

; (Y - 2*Pradius - 2*Pspace) 

WRITE (50, 10) 'The area percentage of particles is', AP 
'X dimension is',X 
'Y dimension is',Y 

'Min spacing between particles is',Pspace 
'Particle radius is', Pradius 
'Particle area is', Parea 
Total number of particles is', Pnum 
'Number of particles along the perimeter 



WRITE(50,20) 
WRITE(50,20) 
WRITE(50,20) 
WRITE(50,20) 
WRITE(50,20) 
WRITE(50, 20) 
WRITE(50,20) 



: is' ,b - 1 

WRITE (50, 20) 'Number of particles in area interior is', a 
1 

WRITE (50, 20) 'Number particles per unit area', NPA 
WRITE (50, 20) 'Average spacing between nearest neighbors 
;is', dpavg 

WRITE(50, 20) 'Min nearest neighbor spacing is', dpmin 
WRITE (50, 20) 'Max nearest neighbor spacing is', dpmax 
PRINT*,' The area percentage of particles is', AP 



PRINT* , 
PRINT* , 
PRINT*, 
PRINT* , 
PRINT* , 
PRINT* , 
PRINT* , 
; 1 

PRINT*, 
PRINT* , 
PRINT* , 
; dpavg 
PRINT* , 
PRINT* , 



X dimension is',X 
Y dimension is',Y 

Min spacing between particles is',Pspace 
Particle radius is', Pradius 
Particle area is', Parea 
Total number of particles is', Pnum 
Number of particles along the perimeter'. 



b - 



Number 

Number 



of particles in area interior', a - 1 
particles per unit area is', NPA 
Average spacing between nearest neighbors is', 

Min nearest neighbor spacing is', dpmin 
Max nearest neighbor spacing is', dpmax 
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10 FORMAT ( //////,T8,A,F10.4) 

20 FORMAT (T8,A,F10.4) 

END 
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APPENDIX C 



PROGRAM RANPART3.FOR 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



PROGRAM RANPART3 



Glossary 

Pi 

Xpos 

Ypos 

XposB 

YposB 

AP 

TPA 

Pspace 

TOTarea 

Parea 

dp 

dp2 

dp2B 

dp3 

dp4 

dist 

dpavg 

X 

Y 

dpmin 

dp2Bmin 

dpmax 

dpinc 

Delta 

Np 

NPA 

Pnum 

R 

Rad(i) 

RadB(i) 

S 

u 

Radmax 

Radmin 

RadBmax 



of variables: 

3 . 14159 

X coordinate of an interior particle 
Y coordinate of an interior particle 
X coordinate of a particle along the border 
Y coordinate of a particle along the border 
Area percentage of the particles 
Total area occupied by the particles 
Min. spacing between particles 
Defined area of particles 
Area occupied by each particle 
Distance between interior particles 
Distance between interior particles 
Distance between an interior particle and a 
border particle 
Nearest neighbor spacing 
Sum of nearest neighbor spacings 
Particle diameter plus min spacing 
Average nearest neighbor spacing 
Horizontal dimension of total area 
Vertical dimension of total area 
Minimum nearest neighbor spacing 
Distance between an interior particle and 
its nearest border particle 
Maximum nearest neighbor spacing 
Size of each range between dpmin and dpmax 
based on the Sturgis Rule 
Minimum value of each range in dpinc 
Number of particles corresponding to each 
Delta 

Number of particles per unit area 

Total number of particles generated 

Generated particle radius from lonormal 

distribution 

Stored particle radius 

Stored border particle radius 

Standard deviation 

Log mean particle radius 

Maximum particle radius 

Minimum particle radius 

Largest border particle 



REAL Pi, Xpos(lOOOO), Ypos(lOOOO), Rad(lOOOO), 
;XposB(300), YposB(300), RadB(300), R, TPA, AP, Pspace, 
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:TOTarea, Parea, dp, TRad, dp2, dp2B, dp2Bmin, u, S, 
:dp3 (10000), dp4, dist, dpavg, X, Y, dpmin, NPA, 

:dpmax, dpinc, g, Delta(50), Np(50) , Pnum, RadBmax, 

: Radmax , Radmin 

* 

INTEGER a, b, c, d, e, f, h, i, j, k, 1, m, n, p, q, nq 

* 

PARAMETER (Pi=3. 14159) 

* 

DATA n, dp4, e, dpmin, dpmax, Np(l), a, b, Parea, TRad, 
; RadBmax, Radmax, Radmin, nr, nq / 

; 1 , 0 . 0 , 0 , 10 . 0 , 0 . 001 , 0 , 1 , 1 , 0 . 0 , 0 . 0 , . 00001 , 

: . 00001 , 1000 , 0 , 0 / 

* 

OPEN (10, file='rp3x.dat' , status='NEW' ) 

OPEN (20,file='rp3y.dat' ,status='NEW'i 
OPEN (30,file='Del3.dat' , status='NEW' ) 

OPEN (40,file='Np3.dat' , status='NEW' ) 

OPEN (50,file='rp3.dat' ,status='NEW' ) 

OPEN (60,file='Rad3.dat' , status=' NEW' ) 

* 

PRINT*,' Enter the area percentage of particles.' 
READ*, AP 

PRINT*,' Enter the length of the X - scale.' 

READ*, X 

PRINT*,' Enter the length of the Y - scale.' 

READ*, Y 

TOTarea = X*Y 

TPA = AP/100*TOTarea 

PRINT*,' Enter the mean particle radius and standard 
tdeviation. . . ' 

READ*, u, S 

PRINT*,' Enter the minimum spacing between particles.' 
READ* , Pspace 

PRINT*,' Enter the total number of particles...' 

READ*, Pnum 

* 

* Determine random position and radius of each particle. 

* 

ISEED = 0 

CALL RNSET (ISEED) 

DO 100 i = 1, Pnum 
175 CONTINUE 
q = 0 

CALL RNLNL (1, U, S, R) 

150 CONTINUE 

q = q + 1 

IF (q .GE. 600) THEN 
nq = nq + 1 
GOTO 175 
ENDIF 



99 



Xpos(i) = RNUNFO *X 
Ypos(i) = RNUNF()*Y 
Rad(i) = R 

* 

* Check for overlap of particles 

* 

DO 200 k = 1, i - 1 

dp = SQRT( (Xpos(i) - Xpos(k))**2 + 

: (Ypos(i) - Ypos(k))**2) 

dist = Rad(i) + Rad(k) + Pspace 
IF (dp .LT. dist) GOTO 150 

200 CONTINUE 

Parea = Parea + Pi*Rad(i)**2 
IF (Parea .GE. TPA) THEN 
Pnum = i 
GOTO 250 
ENDIF 

100 CONTINUE 

250 CONTINUE 

* 

* Segregate particles along the perimeter from those in 

* the interior of the specified area. 

* 

DO 275 d = 1, Pnum 

IF ((Xpos(d) .LE. Rad(d) + Pspace) .OR. 

: (Xpos(d) .GE. X - Rad(d) - Pspace) .OR. 

: (Ypos(d) .LE. Rad(d) + Pspace) .OR. 

; (Ypos(d) .GE. Y - Rad(d) - Pspace)) THEN 

RadBmax = MAX (Rad (d), RadBmax) 

ENDIF 

275 CONTINUE 

DO 280 p = 1, Pnum 

IF ((Xpos(p) .LE. RadBmax + Pspace) .OR. 

; (Xpos(p) .GE. X - RadBmax - Pspace) .OR. 

: (Ypos(p) .LE. RadBmax + Pspace) .OR. 

: (Ypos(p) .GE. Y - RadBmax - Pspace)) THEN 

XposB(b) = Xpos(p) 

YposB(b) = Ypos(p) 

RadB(b) = Rad(p) 
b = b + 1 
ELSE 

Xpos(a) = Xpos(p) 

Ypos(a) = Ypos(p) 

Rad (a) = Rad(p) 

TRad = TRad + Rad (a) 

WRITE(10,*) Xpos(a) 

WRITE(20,*) Ypos(a) 

WRITE(60,*) Rad(a) 
a = a + 1 
ENDIF 

280 CONTINUE 
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* Determine the average distance between nearest 

* neighbors. 

* 

DO 300 j = 1, Pnum - (b - 1) 
dp3(n) = SQRT(X**2 + Y**2) 
dp2Bmin = SQRT(X**2 + Y**2) 

DO 350 c = 1, (b - 1) 

dp2B = SQRT( (Xpos( j) - XposB(c))**2 + 

: (Ypos(j) - YposB(c) ) **2) 

dp2Bmin = MIN(dp2Bmin, dp2B) 

350 CONTINUE 

DO 400 1=1, Pnum - (b - 1) 

dp2 = SQRT( (Xpos (1) - Xpos(j))**2 + 

: (Ypos(l) - Ypos(j))**2) 

IF (dp2 .NE. 0.0) THEN 
dp3(n) = MIN(dp2, dp3 (n) , dp2Bmin) 

ENDIF 

400 CONTINUE 

n = n + 1 

300 CONTINUE 

DO 500 m = 1, Pnum - (b - 1) 
dp4 = dp4 + dp3 (m) 

500 CONTINUE 

dpavg = dp4/(Pnum - (b - 1)) 

* 

* Determine the number of particles falling in each 

* range of nearest neighbor spacings (based on the Sturgis 

* rule) . 

k 

DO 600 f = 1, Pnum - (b - 1) 
dpmin = MIN(dp3(f), dpmin) 
dpmax = MAX(dp3(f), dpmax) 

Radmin = MIN (Rad (f), Radmin) 

Radmax = MAX (Rad (f), Radmax) 

500 CONTINUE 

dpinc = (dpmax - dpmin) /(I + 3 . 3*LOG10 (Pnum - (b - 1) ) ) 

DO 700 g = dpmin, dpmax, dpinc 
e = e + 1 

DO 800 h = 1, Pnum - (b - 1) 

IF ((dp3(h) .GE. g) .AND. (dp3(h) .LT. g + dpinc)) 
: THEN 

Np(e) = Np(e) + 1 
ENDIF 

500 CONTINUE 

Delta(e) = (2*g + dpinc) /2 
WRITE(30,*) Delta (e) 

WRITE (40 ,*) Np (e) / (Pnum - (b - 1)) 

'00 CONTINUE 

NPA = (Pnum - (b - 1))/(X - 2*RadBmax - 2*Pspace) / 

: (Y - 2*RadBmax - 2*Pspace) 
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10 

20 

30 



WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
D) 

WRITE ( 
WRITE ( 
WRITE ( 
1 

WRITE ( 
1 

WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
FORMAT 
FORMAT 
FORMAT 
END 



50, 10) 'The area percentage of particles is', AP 

50.20) 'X dimension is',X 

50.20) 'Y dimension is',Y 

50.20) 'Min spacing between particles is',Pspace 

50.20) 'Mean particle radius is', TRad/ (Pnum - (b 

50.20) 'Total area occupied by particles', Parea 

50.20) 'Total number of particles is', Pnum 

50 , 30) 'Number of particles along the border is',b 

50 , 30) 'Number of particles in area interior is', a 

50. 20) 'Number particles per unit area', NPA 

50.20) 'Average nearest neighbor spacing is', dpavg 

50.20) 'Min nearest neighbor spacing is',dpmin 

50.20) 'Max nearest neighbor spacing is',dpmax 

50.20) 'Max particle radius is',Radmax 

50.20) 'Min particle radius is',Radmin 

,' The area percentage of particles is', AP 

,' X dimension is',X 

,' Y dimension is',Y 

,' Min spacing between particles is',Pspace 

,' Mean particle radius is', TRad/ (Pnum - (b - 1) ) 
,' Total area occupied by particles is', Parea 
,' Total number of particles is', Pnum 
,' Number of particles along the border', b - 1 

, ' Number of particles in area interior' , a - 1 

,' Number particles per unit area is', NPA 
,' Average nearest neighbor spacing is', dpavg 
,' Min nearest neighbor spacing is',dpmin 
,' Max nearest neighbor spacing is',dpmax 
,' Max particle radius is',Radmax 

,' Min particle radius is',Radmin 

(T6,A,T50,F10.4) 

(T6,A,T50,F10.4) 

(T6,A,T46,I9) 
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APPENDIX D. PROGRAM RANPART4.FOR 





PROGRAM 


RANPART4 


* 


Glossary 


of variables: 


* 


Pi 


3.14159 


* 


Xpos 


X coordinate of an interior particle 


* 


Ypos 


Y coordinate of an interior particle 


* 


XposB 


X coordinate of a particle along the border 


* 


YposB 


Y coordinate of a particle along the border 


* 


AP 


Area percentage of the particles 


* 


TPA 


Total area occupied by the particles 


* 


Pspace 


Min. spacing between particles 


* 


TOTarea 


Defined area of particles 


* 


Parea 


Area occupied by each particle 


* 


dp 


Distance between interior particles 


* 


dp2 


Distance between interior particles 


* 


dp2B 


Distance between an interior particle and a 


* 




border particle 


* 


dp 3 


Nearest neighbor spacing 


* 


dp4 


Sum of nearest neighbor spacings 


* 


dist 


Particle diameter plus min spacing 


* 


dpavg 


Average nearest neighbor spacing 


* 


X 


Horizontal dimension of total area 


* 


Y 


Vertical dimension of total area 


* 


dpmin 


Minimum nearest neighbor spacing 


* 


dp2Bitiin 


Distance between an interior particle and 


* 




its nearest border particle 


* 


dpmax 


Maximum nearest neighbor spacing 


* 


dpinc 


Size of each range between dpmin and dpmax 


* 




based on the Sturgis Rule 


* 


Delta 


Minimum value of each range in dpinc 


* 


Np 


Number of particles corresponding to each 


* 




Delta 


* 


NPA 


Number of particles per unit area 


* 


Pnum 


Total number of particles generated 


* 


R 


Generated particle radius from lonormal 


* 




distribution 


* 


Rad ( i ) 


Stored particle radius 


* 


RadB(i) 


Stored border particle radius 


* 


S 


Standard deviation 


★ 


u 


Log mean particle radius 


* 


Radmax 


Maximum particle radius 


* 


Radmin 


Minimum particle radius 


* 


RadBmax 


Largest border particle 


* 


AFl-5 


Area percentage of particles in each of the 


* 




five sub-areas 


* 


Al-5 


Subtotal of particles in each sub-area 
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REAL Pi, Xpos(lOOOO), Ypos(lOOOO), Rad(lOOOO), 
:XposB(300), YposB(300), RadB(300), R, TPA, AP, Pspace, 
:u, TOTarea, Parea, dp, TRad, dp2 , dp2B, dp2Bmin, S, 
:dp3(10000), dp4, dist, dpavg, X, Y, dpmin, NPA, dpmax, 
:dpinc, g, Delta(50), Np(50), Pnum, RadBmax, Radmax, 
tRadmin, Al, AFl, A2, AF2, A3, AF3 , A4 , AF4, A5, AF5 

* 

INTEGER a, b, c, d, e, f, h, i, j, k, 1, m, n, p, q, nq 

* 

PARAMETER (Pi=3.14159) 

* 

DATA n, dp4, e, dpmin, dpmax, Np(l), a, b, Parea, TRad, 
: RadBmax, Radmax, Radmin, Al, A2, A3, A4, A5, nr, nq / 
: 1 , 0 . 0 , 0 , 10 . 0 , 0 . 001 , 0 , 1 , 1 , 0 . 0 , 0 . 0 , . 00001 , 

:. 00001 , 1000 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 , 0 / 

* 

OPEN (10,file='rp4x.dat' , status='NEW' ) 

OPEN (20, f ile='rp4y.dat' , status='NEW' ) 

OPEN (30,file='Del4.dat' , status='NEW' ) 

OPEN (40, file= 'Np4.dat' , status='NEW' ) 

OPEN (50,file='rp4.dat' ,status='NEW') 

OPEN (60,file='Rad4.dat' , status='NEW' ) 

* 

PRINT*,' Enter the area percentage of particles.' 
READ*, AP 

PRINT*,' Enter 5 sub-area percentages of particles' 

READ*, AFl, AF2, AF3 , AF4 , AF5 

AFl = AFl/100 

AF2 = AF2/100 

AF3 = AF3/100 

AF4 = AF4/100 

AF5 = AF5/100 

PRINT*,' Enter the length of the X - scale.' 

READ* , X 

PRINT*,' Enter the length of the Y - scale.' 

READ*, Y 

TOTarea = X*Y 

TPA = AP/100*TOTarea 

PRINT*,' Enter the mean particle radius and standard 
: deviation. . . ' 

READ*, u, S 

PRINT*,' Enter the minimum spacing between particles.' 
READ*, Pspace 

PRINT*,' Enter the total number of particles...' 

READ* , Pnum 

* 

* Determine random position and radius of each particle. 

* 

ISEED = 0 
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175 



150 



* 

* 

* 



200 



CALL RNSET (ISEED) 

DO 100 i = 1, Pnum 
CONTINUE 
q = 0 

CALL RNLNL (1, U, S, R) 

CONTINUE 
q = q + 1 

IF (q .GE. 600) THEN 
nq = nq + 1 
GOTO 175 
END IF 

Xpos(i) = RNUNF()*X 
Ypos(i) = RNUNF()*Y 
Rad(i) = LOG(R) 

Check for overlap of particles 

DO 200 k = 1, i - 1 

dp = SQRT( (Xpos(i) - Xpos(k))**2 + 

: (Ypos(i) - Ypos(k))**2) 

dist = Rad(i) + Rad(k) + Pspace 
IF (dp .LT. dist) GOTO 150 
CONTINUE 

IF ((Ypos(i) .GE. 0.0) .AND. (Ypos(i) .LT. 0.2*Y)) THEN 
IF (A1 .GT. TPA*AF1) GOTO 150 
A1 = A1 + Pi*Rad(i)**2 

ELSEIF ((Ypos(i) .GE. 0.2*Y) .AND. (Ypos(i) .LT. 0.4*Y)) 
:THEN 

IF (A2 .GT. TPA*AF2) GOTO 150 
A2 = A2 + Pi*Rad(i)**2 

ELSEIF ((Ypos(i) .GE. 0.4*Y) .AND. (Ypos(i) .LT. 0.6*Y)) 
:THEN 

IF (A3 .GT. TPA*AF3) GOTO 150 
A3 = A3 + Pi*Rad(i)**2 



ELSEIF ((Ypos(i) .GE. 0.6*Y) .AND. (Ypos(i) .LT. 0.8*Y)) 
: THEN 

IF (A4 .GT. TPA*AF4) GOTO 150 
A4 = A4 + Pi*Rad(i)**2 

ELSEIF ((Ypos(i) .GE. 0.8*Y) .AND. (Ypos(i) .LE. Y) ) 
:THEN 

IF (A5 .GT. TPA*AF5) GOTO 150 
A5 = A5 + Pi*Rad(i)**2 
ENDIF 

Parea = Parea + Pi*Rad(i)**2 
IF (Parea .GE. TPA) THEN 
Pnum = i 



GOTO 250 



ENDIF 



100 CONTINUE 
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CONTINUE 



250 

* 

* Segregate particles along the border from those in 

* the in of the specified area. 

* 

DO 275 d = 1, Pnum 

IF ((Xpos(d) .LE. Rad(d) + Pspace) .OR. 

: (Xpos(d) .GE. X - Rad(d) - Pspace) .OR. 

: (Ypos(d) .LE. Rad(d) + Pspace) .OR. 

: (Ypos(d) .GE. Y - Rad(d) - Pspace)) THEN 

RadBmax = MAX (Rad (d) , RadBmax) 

ENDIF 

275 CONTINUE 

DO 280 p = 1, Pnum 

IF ((Xpos(p) .LE. RadBmax + Pspace) .OR. 

: (Xpos(p) .GE. X - RadBmax - Pspace) .OR. 

: (Ypos(p) .LE. RadBmax + Pspace) .OR. 

: (Ypos(p) .GE. Y - RadBmax - Pspace)) THEN 

XposB(b) = Xpos(p) 

YposB(b) = Ypos(p) 

RadB(b) = Rad(p) 
b = b + 1 
ELSE 

Xpos(a) = Xpos(p) 

Ypos(a) = Ypos(p) 

Rad (a) = Rad(p) 

TRad = TRad + Rad (a) 

WRITE(10,*) Xpos(a) 

WRITE(20,*) Ypos(a) 

WRITE(60,*) Rad(a) 
a = a + 1 
ENDIF 

280 CONTINUE 

* 

* Determine the average distance between nearest 

DO 300 j = 1, Pnum - (b - 1) 
dp3(n) = SQRT(X**2 + Y**2) 
dp2Bmin = SQRT(X**2 + Y**2) 

DO 350 c = 1, (b - 1) 

dp2B = SQRT( (Xpos( j) - XposB(c))**2 + 

: (Ypos(j) - YposB(c) ) **2) 

dp2Bmin = MIN(dp2Bmin, dp2B) 

350 CONTINUE 

DO 400 1=1, Pnum - (b - 1) 

dp2 = SQRT( (Xpos(l) - Xpos(j))**2 + 

: (Ypos(l) - Ypos(j))**2) 

IF (dp2 .NE. 0.0) THEN 
dp3 (n) = MIN(dp2, dp3(n), dp2Bmin) 

ENDIF 

400 CONTINUE 
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n = n + 1 

300 CONTINUE 

DO 500 m = 1, Pnum - (b - 1) 
dp4 = dp4 + dp3 (m) 

500 CONTINUE 

dpavg = dp4/(Pnum - (b - 1)) 

* 

* Determine the number of particles falling in each range 

* of nearest neighbor spacings (based on the Sturgis 

* rule) . 

* 

DO 600 f = 1, Pnum - (b - 1) 
dpmin = MIN(dp3(f), dpmin) 
dpmax = MAX(dp3(f), dpmax) 

Radmin = MIN(Rad(f), Radmin) 

Radmax = MAX (Rad (f) , Radmax) 

600 CONTINUE 

dpinc = (dpmax - dpmin) /(I + 3 . 3*LOG10 (Pnum - (b - 1) ) ) 
DO 700 g = dpmin, dpmax, dpinc 
e = e + 1 

DO 800 h = 1, Pnum - (b - 1) 

IF ((dp3(h) .GE. g) .AND. (dp3 (h) .LT. g + dpinc)) 

: THEN 

Np(e) = Np(e) + 1 
ENDIF 

800 CONTINUE 

Delta(e) = (2*g + dpinc) /2 
WRITE (30,*) Delta (e) 

WRITE(40,*)Np(e)/(Pnum - (b - 1) ) 

700 CONTINUE 

NPA = (Pnum - (b - 1))/(X - 2*RadBmax - 2*Pspace)/ 

: (Y - 2*RadBmax - 2*Pspace) 

WRITE (50, 10) 'The area percentage of particles is', AP 
WRITE (50, 50) 'The sub-area percentages of particles 
: are: ' 

WRITE (50, 40) AF1,AF2,AF3,AF4,AF5 
WRITE(50, 20) 'X dimension is',X 
WRITE(50,20) 'Y dimension is',Y 

WRITE(50, 20) 'Min spacing between particles is',Pspace 
WRITE (50, 20) 'Mean particle radius is', TRad/ (Pnum - (b 
1)) 

WRITE(50, 20) 'Total area occupied by particles', Parea 
WRITE(50, 20) 'Total number of particles is', Pnum 
WRITE (50, 30) 'Number of particles along the border is',b 
1 

WRITE (50, 30) 'Number of particles in area interior is', a 
1 

WRITE (50, 20) 'Number particles per unit area', NPA 
WRITE(50,20) 'Average nearest neighbor spacing is', dpavg 
WRITE(50, 20) 'Min nearest neighbor spacing is', dpmin 
WRITE (50, 20) 'Max nearest neighbor spacing is', dpmax 
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10 

20 

30 

40 

50 



WRITE(50,20) 'Max particle radius is',Radmax 
WRITE(50, 20) 'Min particle radius is',Radmin 
PRINT*,' The area percentage of particles is', AP 
PRINT' (T6 , A, T45 , 5F4 . 1) ' , 'The sub-area percentages of 
: particles are ' , AFl , AF2 , AF3 , AF4 , AF5 

PRINT' (T6,A,T45,5F6.2) ' , 'The area occupied by particles 
: in each sub-area is...',Al, A2, A3, A4, A5 



PRINT*,' X dimension is',X 

PRINT*,' Y dimension is',Y 

PRINT*,' Min spacing between particles is',Pspace 

PRINT*,' Mean particle radius is', TRad/(Pnum - (b - 1) ) 

PRINT*,' Total area occupied by particles is', Parea 
PRINT*,' Total number of particles is', Pnum 
PRINT*,' Number of particles along the border', b - 1 

PRINT*,' Number of particles in area interior', a - 1 

PRINT*,' Number particles per unit area is', NPA 
PRINT*,' Average nearest neighbor spacing is', dpavg 
PRINT*,' Min nearest neighbor spacing is',dpmin 
PRINT*,' Max nearest neighbor spacing is',dpmax 
PRINT*,' Max particle radius is',Radmax 

PRINT*,' Min particle radius is',Radmin 

FORMAT(T6,A,T50,F10.4) 

FORMAT(T6,A,T50,F10.4) 

FORMAT (T6, A, T45, 110) 

FORMAT(5F6.2) 

FORMAT (A) 

END 
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